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Abstract 

The proton-proton {pp) fusion cross-section found at the heart of solar models is unconstrained 
experimentally and relies solely on theoretical calculations. Effective field theory provides an 
opportunity to constrain the pp cross-section experimentally, however, this method is complicated 
by the appearance of two-nucleon effects in the form of an unknown parameter Li^a- We present a 
method to constrain Li^a using the Standard Solar Model and helioseismology. Using this method, 
we determine a value of Li^a = 7.0 fm^ with a range of 1.6 to 12.4 fm^. These results are consistent 
with theoretical estimates of Li^a ~ 6 fm'^. 

PACS numbers: 21.45. +v,23.40.Bw,26.20.+f,96.60.Ly 
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I. INTRODUCTION 



The Standard Solar Model (SSM) is a model of the sun constructed by numerically solving 
the equations of stellar structure and evolution constrained by the observed composition, 
luminosity, and radius of the sun at the current age. This model yields acoustic oscillation 
frequencies (p-modes) that match the observed p-mode oscillation spectrum of the sun to 
better than 0.3% [|1|. The uncertainties associated with the observed oscillation frequencies 
are an order of magnitude smaller. The excellent agreement between the model and observed 
frequency spectra implies that the run of sound speed predicted by the model is close to 
that of the sun which, in turn, suggests that the physics used to construct the model of the 
sun is accurate. 

The individual p-mode frequencies extend to distinct depths in the sun, and hence can 
be used to study the physics of specific regions. Of interest to the work described here 
are the low-/ (where / corresponds, in spherical harmonic nomenclature, to the number of 
nodes in azimuth) p-mode frequencies that penetrate deep into the core of the sun. These 
modes, grouped in combinations that cancel out surface dependencies, can be used to test 
the important physics of the core. One of the greatest sources of uncertainty in the physics 
of the solar model core is the pp-fusion cross-section. The reaction is unconstrained by 
laboratory-based experiments, relying completely on precise calculations from theoretical 
nuclear physics. Owing to the precise agreement between the observed and predicted os- 
cillation frequencies that already exists, it is possible to use the oscillation frequencies to 
constrain the pp cross-section, that is, to determine the range of cross-sections that yield 
oscillation frequencies within the known uncertainties of the solar model physics and con- 
straints. To date, this sort of constraint has been imposed only through simplifications of 
the standard solar model 0, and here we use a fully-developed stellar evolution code, to be 
discussed later. 

The cross-section for pp — > de'^i'e has been studied extensively P, ||, 0, H, |Tl 



T2| , |T^, [T^ and contains two components. First, an 'Impulse Approximation' (lA) contribu- 



tion where the weak interaction takes place on a single proton, as shown in Fig. |I]a). This 
comprises more than 95% of the cross-section at energies of solar interest and is very well 
understood because it concerns only the weak properties of a proton and well-known proton- 
proton scattering physics. The second component, and the remainder of the cross-section. 
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FIG. 1: The process pp de'^Ue- o-) The Impulse Approximation (lA), where the weak interaction 
takes place on one proton while the other proton is just a spectator. The final state, represented 
by the grey bubble, is the resulting neutron and proton being bound as a deuteron. b) A two- 
nucleon process where the protons are undergoing a strong interaction at the same time as the weak 
interaction (represented by the black circle). 

is unconstrained by experiment and has led theorists to try to improve our understanding 
of this process. Conventional potential model calculations will introduce so-called meson- 
exchange currents; essentially two-nucleon effects where the weak interaction occurs while 
the two nucleons are interacting with each other. An example of such a process is shown in 
Fig. |l|b). These two-nucleon effects introduce physics that cannot be constrained directly by 
pp elastic scattering experiments, and the parameters of this extended physics must be con- 
strained using other methods. Few other methods exist. Neutrino-deuteron scattering could, 
in principle, constrain the two-nucleon effects. Experiments involving neutrino-deuteron 
scattering are, however, quite difficult. Reactor antineutrino-deuteron experiments can be 
performed, where the experimental uncertainties are of order 10-20%. Nonetheless, they 
can be used to infer a constraint on the two-nucleon matrix elements |T^. One interesting 
method is to use tritium /3-decay to constrain the two-nucleon matrix elements needed in 
pp —>■ de^Vf. [r3|. The tritium half-life is well known, and the same two-nucleon effects occur 



here as in pp fusion. However, there are three nucleons in tritium, hence we cannot be cer- 
tain that this more complicated structure does not, in turn, complicate the weak interaction 
physics further. 

If we were to rely purely on theory to predict the two-nucleon component of the pp 
fusion cross-section, we still have one more problem. How do we define the uncertainty 
in a theoretical calculation? In quantum field theories, such as quantum electrodynamics, 
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there is a perturbative expansion parameter that dictates the size of contributions at each 
order in perturbation theory and thus the size of any neglected remainder. The remainder 
represents the uncertainty in the theory. However, conventional nuclear physics calculations 
use Schrodinger's equation and nucleon-nucleon potentials to calculate initial and final-state 
wavefunctions and, in turn, calculate the matrix elements of relevant transition operators. 
In this approach, no scheme exists for estimating the size of contributions from unknown or 
neglected physics. An alternative to using wavefunctions to study our problem is to use an 
effective field theory (EFT). An EFT has many benefits, but there are drawbacks. If the EFT 
is perturbative (i.e., there is a small expansion parameter), then one can compute to a specific 
order in the parameter and then make direct statements of the scale of theoretical uncertainty 
of the problem as determined by the relative size of the next term in the expansion - just 
as one would do in quantum electrodynamics. This is particularly important given that 
EFTs typically have an infinite number of terms, and the expansion parameter is required 
to ensure that only a finite number of terms are needed for any given calculation. In the 
particular case of nuclear EFTs, specifically ones where all interactions are of zero-range (i.e., 
no pions), the expansion parameter is denoted by e = Q/M^r, where Q represents a generic 
momentum in the problem and M^^ is the pion mass. An added benefit of this approach 
is that field theories pose no problems for conservation laws such as electromagnetic gauge 
invariance, which are often difficult to preserve when wavefunctions are used. However, 
EFTs cannot escape the problem that all parameters in an EFT must be constrained by 
experiment, including the various weak interaction parameters. 



II. PROTON-PROTON FUSION 



There have been a number of works dealing with pp-fusion in EFT [|16], |T^, [T8|, using 



the EFT without pions developed in refs. |]T9|, ^ |2^, As with the traditional 

calculations, the strength of the weak interaction on an interacting pair of nucleons is un- 



known and this ignorance can be specified through one parameter, which we label Li a 



Li A represents the strength of the interaction shown in Fig. Pd). It also appears in the 
various z/(z7) — d breakup reactions and a precise measurement of one of those reactions 
would also constrain pp-fusion. However, precise measurements do not exist (but have been 



proposed [p6| ). 
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It is also interesting to consider the inverse of this argument. A constraint on pp-fusion 
would lead to a constraint on z/(z7) — d breakup reactions. The recent results from the 
Sudbury Neutrino Observatory (SNO) |2^, |28[ rely on the reaction v^d — > e~pp to detect 
and extract the solar neutrino flux, hence the theoretical uncertainties in the cross-section 
certainly affect their results. In spite of the importance of SNO, not as much attention has 
been paid to these breakup reactions, particularly the pair 

fed -> e'pp , (1) 
u^d v^pn . (2) 

The first reaction commonly known as the Charged- Current (CC) reaction, is sensitive to 
the electron neutrino (z/g) flux only. The second Neutral- Current (NC) reaction is equally 
sensitive to electron, muon, and tau neutrinos (denoted by Ux), and provides the key test 
for neutrino oscillations. The CC reaction is very closely related to pp-fusion, but the finite 
energy /momentum transfers in the CC and NC reactions make conventional calculations 
much more difficult than the (relatively) simpler pp-fusion calculation. However, both pp- 
fusion and deuteron breakup can be calculated using the same EFT in a straightforward 
manner and the connection between the two is discussed elsewhere |T^. Here, we wish to 
explore the possibility of using the high-precision associated with helioseismology to impose 
some additional constraints on Li ^. In turn, the implications of these constraints can be 
seen for v — d breakup and for SNO. 

Cross-sections of astrophysical interest are written as 

a(E) = ^exp(-27rr?), (3) 

separating the nuclear physics {S{E)) from the longer range Coulomb effects (the Gamow 
factor exp(— 27r?7)). In the particular case of pp-fusion, we are interested in Spp{0). Adel- 
berger et al. parameterize 5*^^(0) as 



S„(0) = 4.00 X lT-^M.VbO£^^Y(J^\('iil^S^^M— 

' \ 3073 s / V6.92/ V1.2654/ VO. 144/ V 1.01 



(4) 



where the parameters are described in ref. [2£]. Of importance here are A and 5. Their 



parameterization uses an explicit separation of lA and two-nucleon effects into A and 5 
respectively. In this separation the total matrix element for the fusion process is given by 
A(l + 5). This separation is artificial and is not appropriate when a fully self-consistent 



nuclear physics calculation is performed. Consequently, we set 5 = and incorporate all 
effects directly into A. 

In EFT, it was found that [g 



A = 2.58 + 0.011 



I^) - 0.0003 (^). (5) 



where Li^a parameterizes two-nucleon physics at next-to-leading order (NLO), and Ki a 
parameterizes two-nucleon physics at next-to-next-to-next-to- leading order (NNNLO). No 
new parameters appear at the intermediate order (NNLO). Naively, each order of the ex- 
pansion will be (naively) 30% of the previous order. So, a calculation to N^LO should have 
a precision of 1%. In practice, low-energy calculations of this type are usually more rapidly 
convergent. Li^a and Ki a are both unknown, and need to be constrained experimentally. 



Dimensional analysis as developed in refs. [24, 25 1 favors 
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I^i.aI ~ 72 ^ 6 fm% 

[m^ - 7) 

\Kia\ ~ —TT^ ^3 ~ 20 fm^ . (6) 

mi [m^ - 7) 

It can be shown that Ki^a would contribute at an order much less than 1%, so we neglect 



it 18 . 



III. HELIOSEISMOLOGY AND THE STANDARD SOLAR MODEL 

The solar models for this investigation were constructed using the Yale stellar evolution 
code (YREC) in its non-rotating configuration |^, ^ |31|. The code implements the most 



current constitutive physics available. The equation of state tables are provided by the 



OPAL group |32|, The opacities in the interior are interpolated from the OPAL opacity 



tables P3], and the lower temperature opacities, near the surface, are interpolated from the 



Alexander and Ferguson tables. 



The nuclear reaction cross-sections have recently been updated (compared to |[T|) to 
correspond to the latest agreed upon values set at the 1998 Seattle Workshop on Solar 
Fusion Reactions. The nuclear energy generation routines are nearly identical to those used 



by Bahcall et al. |^ who, themselves, use a variation of YREC for their solar model 



calculations. To facilitate quick adjustments of S-factors for each reaction, YREC uses 



Bahcall |^ values multiplied by a scaling factor 5'(0)/5'*(0), where 5(0) is the revised 



S-factor, and 5'*(0) is the original Bahcall value [0. For the results presented in this 
investigation, we adjust the scaling factor corresponding to the pp —>■ de^Ve reaction. The 
models were evolved from the zero age main sequence (ZAMS, the point in the evolution of 
the star where the nuclear power first dominates gravitational power) in 100 equally spaced 
time steps to 4.55 Gyr. The meteoritic age of the Sun, that is the age implied from the age 
of the oldest meteorites, is 4.52±0.04 Gyr [^]. Here we use 4.55 Gyr for the age of the solar 
model, which lies within the uncertainty range of the meteoritic age, to facilitate testing and 
comparison of our models to solar models that we and others have calculated in the past. 
The solar radius, Rq = 6.958 x 10^" cm, and the solar luminosity, Lq = 3.8515 x 10^^ erg s~^, 
are known to better than ±0.05%. The mass of the Sun is 1.989 x 10'^^ g with a relative 
uncertainty of ±0.02% [R^. The mixing length parameter, a free parameter in the mixing 



length theory used to model convective energy transport, and the helium abundance, 
an element whose abundance cannot be determined from spectroscopic measurements of the 
photosphere of the Sun, are adjusted to produce models that match each other and the 
above-stated values of Rq and Lq to one part in 10^. 

In order for the surface composition of the (evolved) solar model to match the observa- 
tionally determined abundances, we set the initial (ZAMS) homogeneous composition for 
hydrogen, helium and metals (A>4) by mass fraction to 0.70633, 0.27367, and 0.02000, re- 
spectively. The effects of gravitational settling of helium and metals, which alters the surface 
abundances as the models evolve, are modeled using the diffusion formulation described by 



Bahcall et al. [^. Other modeling details are as described in [|l|. 

To constrain the range for Li^a, a relationship between Li^a and Spp is established. Ne- 
glecting the higher-order effects of Ki a in eq. |^, 

A = 2.58 ± 0.011 (^) (7) 



and applying the appropriate values from to eq. ^ gives 



Spp{0) = 0.56663 X 10-^^(A^)MeVb . (8) 

The combination of eqs. and § establishes a relationship between Li^a and Spp{0), given 
by 

Spp{0) = 3.7717 X 10-^^ (l + 8.53 x 10~^ {j^) + ^'^^ ^ ^^'^ (t^)1 MeVh . (9) 
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FIG. 2: Five curves representing p-mode frequencies at 1000 fiHz, 2000 ^Hz,..., 5000 fiHz in a plot 
of the radius fraction of the inner turning point versus azimuthal order. The inner turning point 
approaches the core as I decreases. 

Using the above relationship and the current theoretical value for Spp{0), the best determined 
value for Li a may be calculated. By varying the Spp{0) value within YREC in accordance 
with values of Li^a, solar models are produced corresponding to each variation. 

The resultant models are then read into a nonradial, nonadiabatic stellar oscillation pro- 
gram where the oscillation frequencies are calculated. We adopt the usual nomenclature 



describing different orders of acoustic oscillations, or p-modes, with the radial order of the 
mode given by n and the azimuthal order of the mode given by /. Only the lowest /-valued 
p-modes (i.e., / = 0, 1) penetrate into the central nuclear burning regions of the sun. The ra- 
dius to which a mode penetrates, the inner turning point r^, is related to the cyclic frequency 



(jj and the sound speed c{rt) by [^| 



^ V^-'-' . (10) 

UJ 

This behaviour is shown in Fig. ^ 

The physics used to model the structure near the surface of the solar model is more 
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uncertain than the physics of the core. Indeed, comparisons of j9-mode frequencies that 
probe the surface reveal greater discrepancies between model and observations than in the 
core. It is believed that the mixing length approximation used to model convection is the 
principal modeling uncertainty. The frequencies of all observed p-modes are affected by this 
region. It is possible, though, to select pairs of p-modes that have identical eigenfunction 
shape near the surface, and distinct shape in the deep interior, and use them in combination 
to cancel out the surface layer effects. This is achieved via the small frequency spacing, 
5z/(n, /), which is defined as the frequency difference between an z/(n, /) and an z/(n — 1, / + 2) 
p-mode. For / = and 1, the 5i/(n, /) is well known as one of the most sensitive helioseismic 
diagnostics of the central-most regions of the Sun. 



IV. RESULTS 



The best value for Li^a is determined by using the current value for Spp{0). From 
the current value for 5*^^(0) is 

Spp{0) = 4.00(1 ± 0.007+°:°20) ^ io-25MeVb (11) 



Using eqs. ^and |TT|, the determined value for the unknown counterterm is found to be, Li a 
= 7.0 fm^. 

To determine the range in Li a using helioseismology, plots of (5z/(n, /) versus adiabatic 
oscillation frequencies are created. To enhance the analysis, Si>{n, I) differences are used 
instead of 5z/(n, /). These differences are formed by taking the 5v{n,l) results from the 
individual models and subtracting the observed 5z/(n, /) for the Sun; Sumodeiin, I) - Suq^ti, I), 
where Si>Q{n,l) values are calculated from observed solar oscillation frequencies. 

Low-/ small spacings have been published for a number of different high quality solar 
p-mode observations, including GONG (Global Oscillation Network Group; [^, and 
BiSON (Birmingham Solar Oscillation Network; ^^). Here we choose the BiSON data set 
because the uncertainties are low and because we have used this data set before and are 
familiar with it. 

A plot of 5z/(n, I) differences versus adiabatic frequencies for the reference SSM is shown 
in Fig. 1^. In Fig. ^, the thick solid line represents the reference SSM model result. The error 



bars correspond to the uncertainties stated for the observed BiSON results for the Sun [44 
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FIG. 3: The small spacing frequencies 5v{n,l) for the reference SSM, shown relative to the ob- 
served values from BiSON. The error bars represent the uncertainties in the BiSON data. The 
discrepancies between the model and observations are largest at high frequencies, which are most 
sensitive to surface effects. The reference SSM uses a value for Li^a of 7.0 fnt", and is shown by 
a solid black line. The shaded area represents the variation of Li^a between 1.6 and 12.4 fi^ ■ 



The low-/ small spacing frequencies of our reference SSM lie almost entirely within the 
uncertainty range of the observationally derived small spacings. This suggests that the 
interior structure of our model is a very close fit to the actual sun. Which, in turn suggests 
that the physics specific to this region are accurate, including the age estimate, the chemical 
abundance and diffusion processes, the opacities, the equation of state, and the nuclear 
reaction rates. That is, no modifications to the current SSM are required to make it agree 
with the most accurate p-mode data we currently have available to us. It is comforting that 
the simplest, non-ad-hoc, physics has produced an accurate match to the observed ]9-mode 
data. 

The analysis described above was repeated for each, nonstandard, model corresponding 
to a variation of Li a with respect to the reference SSM. A quantitative estimate of the range 
of Li A is established by comparing the distribution of the frequency results for various Li^a 
with respect to the observed BiSON error bars. To extract this quantitative result, the 
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frequency range of ~2500 to ~3000 /iHz was studied since this frequency range for the SSM 
showed an excellent fit to the observed BiSON data. 

For each model results, average (5z/(n, /) difference values were calculated within the fre- 
quency range of ~2500 to ~3000 /iHz. The average observed error within this frequency 
range was calculated to be ±0.066 /xHz. Average 5u{n, I) differences are plotted versus the 
corresponding Li^a value for both the upper and lower range. These plots show a linear re- 
lation between the average 5z/(n, /) differences and Li^a- From this relation and the average 
observed error values, the values of Li ^ corresponding to the upper and lower limits of the 
error bars are calculated. The upper limit is determined to be 12.4 fm^. The lower limit is 
determined to be 1.6 fm^. The range of allowed values is shown as the shaded area in Fig. |^. 
Therefore the value and effective range for Li a is determined to be 

Li,A = 7.0 ± 5.4 fm^ (12) 

A further result from the determined SSM is the calculation of the theoretical total 
neutrino flux of 4.93 x 10^ cm^^s^^. This result is consistent with the total observed 
neutrino flux from SNO of (5.09lJ]i^(stat)+!|:43(syst)) x 10^ cm-^g-i |2|. 

We have tacitly assumed that the uncertainty in the model-computed p-mode frequencies 
is comparable to, or less than, the error bars associated with the observed frequencies. This 
is not too unreasonable since the model-computed frequencies do lie within the uncertainty 
error bars of the observations. In order to more formally justify this assumption, we must 
carry out an extensive uncertainty analysis of the solar model physics and their impact on 
the computed p-mode frequencies. This important task, upon which all similar helioseismic 
tests depend, is in progress and will be reported at a later time. 



V. CONCLUSIONS 



Using the relationship between effective field theory, nuclear cross-sections and the accu- 
racy of helioseismology, the unknown counterterm, Li^Ai is determined to have a value of 7.0 
fm'^ with a range of 1.6 to 12.4 fm^. This result for Li^a is consistent with the theoretical 
value determined by Butler & Chen |^ using dimensional analysis of Li^a ~ 6 fm^. The 
result and range is also consistent with other evaluations of Li a which include 5.6±2.0 fm^ 
and 6.5±2.4 fm^ p . 
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The standard solar model determined for this work shows remarkably excellent agreement 
with the latest BiSON oscillation frequency observations. Theoretical total fluxes for 
neutrinos from the standard solar model result also show excellent agreement with total 
observed fluxes reported from the Sudbury Neutrino Observatory. 
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